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Abstract 

By generating function based on the Jackson's q-exponential function 
and standard exponential function, we introduce a new q-analogue of Her- 
mite and Kampe-de Feriet polynomials. In contrast to standard Hermite 
polynomials, with triple recurrence relation, our polynomials satisfy mul- 
tiple term recurrence relation, derived by the q-logarithmic function. It 
allow us to introduce the q-Heat equation with standard time evolution 
and the q-deformed space derivative. We found solution of this equation in 
terms of q-Kampe-de Feriet polynomials with arbitrary number of moving 
zeros, and solved the initial value problem in operator form. By q-analog 
of the Cole-Hopf transformation we find a new q-deformed Burgers type 
nonlinear equation with cubic nonlinearity. Regular everywhere single and 
multiple q-Shock soliton solutions and their time evolution are studied. 
A novel, self-similarity property of these q-shock solitons is found. The 
results are extended to the time dependent q-Schrodinger equation and 
the q-Madelung fluid type representation is derived. 



1 Introduction 

It is well known that the Burgers' equation in one dimension can be reduced 
via the Cole-Hopf transformation to the linear heat equation. It allows one to 
solve the initial value problem for the Burgers equation and get exact solutions 
in the form of shock solitons and their scattering. In the present paper we intro- 
duce the differcntial-q-differcnce Burgers type equation with cubic nonlinearity 
which includes the standard time derivative and the q-deformed space deriva- 
tive. By using the q-Cole-Hopf transformation, this nonlinear equation can be 
linearized in terms of the q-Heat equation with standard time evolution and 
q-different space derivative. Based on the Jackson's q-exponential function and 
the standard exponential function we introduce a new q-analog of Hermite and 
Kampe-de Feriet polynomials, representing moving poles solution for the q- 
Burgers equation. Then we derive the operator solution of the initial value 
problem (I VP) for the q-Burgers equation in terms of the IVP for the q-heat 
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equation. We construct several particular solutions of our q-Burgers type equa- 
tion in the form of regular everywhere q-shock solitons with regular time evo- 
lution free of singularities. It turns out that the static q-shock soliton solution 
of our equation shows remarkable self-similarity property in space coordinate x. 
By extending our results to the complex domain we introduce the time depen- 
dent Schrodinger equation with q-deformed dispersion and the complex wave 
function. As a solution of this equation we get the set of complex q-Kampe- 
de Feriet polynomials. By the complex q-Cole-Hopf transformation we obtain 
the complex q-Burgcrs-Madelung equation as a coupled two fluid system with 
complex velocity function. 



2 q-Hermite Polynomials 

We define a q-analog of Hermite polynomials by the generating function with 
the Jackson's q-exponential function |3j and standard exponential function as 

2 ^ t N 

e~* e q ([2] g tx) = $J ^M^uTTr W 
where the Jackson's q-exponential function is defined by 
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[l] q [2]q...[n] g and q-number pQ, 
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From the defining identity ((T|) it is not difficult to derive for the q-Hermite 
polynomials an explicit sum formula 

HN{x ' q) - ^ k\[N-2k] q \ • (2) 

This explicit sum makes it transparent in which way our polynomials Hff{x\q) 
q-extended the Hn(x) and how they are different from the known ones in lit- 
erature. By ^-differentiating the generating function (|TJ) with respect to x we 
derive the two-terms recurrence relation 

D x H N (x;q) = [2] q [N} q H N _ 1 [x;q), (3) 

where the q-derivative is defined as [1] 

D,m = nf-H*\ ( 4) 

[q - l)x 
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And by standard differentiating the generating function (TTJ) with respect to 
t and using the next evident equality 

n=0 L lq 

we obtain another two-term recurrence relation 

A 

dx 



(x— - N)H N (x; q) = 2[N] q [N - l] q H N _ 2 (x; q). (5) 



By standard differentiating the generating function ([T]) with respect to t and 
using definition of the q-logarithmic function [5] 

wi + ^) = E ( L ■ 

where q > 1,0 < \z\ < q and the property 



d ( az \ Ln q (l — az) 

—Ln e„ 1 1 — 



dz \1 — q J (q—l)z 

we derive the N-term recurrence relation formula 



"\ ! =^A^f{[2} q xH N (x-,q)-2[N} q H N ^(x-,q)- 

[A:] g ![iV-fc + l] g 



(« - l)^]^^;?) + [2} q [N] q \ £ ( rw.rv-t.uii > 



When g — > 1 this multiple term recurrence relation for q-Hermite polynomials 
reduces to the three-term recurrence relation for the standard Hcrmite polyno- 
mials 

H N+1 (x) = 2xH N {x) - 2NH N - 1 {x). 

Substituting ([3]) into N-term recurrence relation formula we get the operator 
representation 

H N+1 (x;q)= (7) 

l^Jk ( [21sI _ , JL + ( , _ ^ )D , + 1 y^n;) h nM 

By the recursion, starting from n — and Hq(x; q) — 1 we have next represen- 
tation for the q-Hermite polynomials 

H N+1 (x;q) = (8) 
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In the limit q — > 1 case this product formula is reduced to the known one 



H N {x) = (2x - ^-) N ■ 1 

We note that the generating function and the form of our q-Hcrmite poly- 
nomials are different from the known ones in the literature, [2], [5], [4], [6], 
[9]. Moreover, instead of three-term recurrence relation we have multiple term 
recurrence relation, which shows that our q-Hermite polynomials are different 
from the known ones for orthogonal polynomial sets [7]. 

The first few q-Hcrmite polynomials are 

H (x;q) = 1, 
Hi{x;q) = [2] q x, 



H 2 (x;q) = [2] 2 q x 2 - [2]„ 



H 3 (x;q) = [2] q x 6 -[2} 2 [3] q x, 

H i (x;q) = [2} q x i -[2} 2 q [3} q [4} q x 2 + l[A} q l. 
When q —} 1 these polynomials reduce to the standard Hermite polynomials. 

2.1 q-Difference Equation 

Applying D x to both sides of ^ and using recurrence formula (|5|) we get the 
q-difference-differential equation for the q-Hermite polynomials 

2D 2 q H N (x; q) - [2} 2 q x—H N (x; q) + [2} 2 q NH N (x; q) - 0. (9) 

In q — > 1 limit it reduces to the second order linear differential equation for 
the standard Hermite polynomials 



-—rH N (x) - 2x— 

ax dx 



2 H N (x) - 2x—H N {x) + NH N (x) = 0. 



3 q-Kampe-de Feriet Polynomials 

We define the q-Kampe-de Feriet polynomials as 

H N (x,ut;q) = {-ut)^H N ^J= ,<l] ■ (KM 
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so that from N-term recurrence relation for q-Hcrmitc polynomials we obtain 
N-term recurrence relation formula for q-Kampe-dc Fcrict polynomials 

H N+1 (x,vt;q) = ^ffi [xH N (x, vt; q) + 2vt[N] q H N _ 1 (x, vt; q) 

This can also be written in the operator form as 

H N+ i (x, vt; q) = 
[N + l] q 

x -r t -rr, — ;u x t 

1=2 



N + 1 



N 



[2], 



2 )^ + E' 



(ii) 

H N (x,vt;q) 



By the recursion, starting from n — and Hq(x, vt; q) = 1 we have next operator 
representation for the q-Kampe-de Feriet polynomials 



H N+1 (x,vt;q) = 

N 



(12) 



TT [fc + lk 

In g — > 1 case we have 



1-5 



™ n _ „2\fe T fe+i 



[2], 



fe=2 



# w (x) = (x + 2vt— ) N ■ 1. 

eta 

Then the first few q-Kampe-de Feriet polynomials are 

H (x,vt;q) = 1, 

-H"i(a;, i4; g) = x, 



i? 2 (a;, ^; 9) = x 2 + [2] q vt, 



H 3 (x, vt; q) — x 3 + [2} q [3} q vtx, 



H 4 (x,vt;q) =x i + [3] q [4] q vti 



[4] 



i- ,,2+2 



vH z . 



When these polynomials reduce to the standard Kampe-de Feriet poly- 

nomials. 
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4 q-Heat Equation 

We introduce the q-Heat equation 

(d t - vDl)<t>(x,t) = (13) 

with partial q-derivative with respect to x and with partial standard derivative 
in time t. One can easily see that 

4> k (x,t) = e vhH e q [hx) 



is a plane wave solution of (|13[) . By expanding this in terms of parameter k 

[N] 



2 ^ k N 

(x,t) = e vh ^(kx) = ]T Hufavtiq)-^ (14) 

JV=0 



we get the set of q-Kampe-de Feriet polynomial solutions for the q-Heat equation 
(fl~3"|) . From the defining identity (fT4")l is not difficult to derive an explicit sum 
formula for the q-Kampe de Feriet polynomials 

H N (x,vt; q )=^2 kl[N _ 2k]ql ■ ( 15 ) 

4.1 Operator Representation 
Proposition 1 

e~^ Dx e q {[2] q xt) = e~ t2 e q ({2} q xt). (16) 
Proof 1 By q- differentiating the q- exponential function with respect to x 

D n x e q {[2] q xt) = ([2}t) n e q ([2] q xt), (17) 
and combining then to the sum 

^ a ~Dl n e q ([2] q xt) = £ [[ lq \ ° e q ([2] q xt), (18) 

n=0 U - n=0 n - 

we have relation 

e aD * e q ([2] q xt) = e a ^ 2 e q ([2] q xt). (19) 



By choosing a = — l/[2] q we get the result 1116] 
Proposition 2 



H N (x;q) = [2]fe ™t Dx x N . (20) 



Proof 2 The right hand side of M6)) is the generating function for the q-Hermite 
polynomials {I]). Hence, equating the coefficients of t n on both sides gives the 
result. 
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Corrollary 1 If function f(x) is expandable to the formal power series f(x) = 
Sat=o a NX N then we have next q-Hermite series 

~W D *f( \ H N (x;q) 
e n " f(x) = }^ a N [9lAf ■ (21) 

N=0 [ J 9 

5 Evolution Operator 

Following similar calculations as in Proposition I we have the next relation 

e vtD *e q {kx) = e vtk2 e q {kx). (22) 

The right hand side of this expression is the plane wave type solution of the q- 
heat equation (1131) . Expanding both sides in power series in fc and equating the 
coefficients of k N on both sides we get q-Kampe de Feriet polynomial solutions 
of this equation 

H N (x,ut;q) =e vtD *x N (23) 

as solution of the I.V.P. with <j>(x, 0) = x N . 

Consider an arbitrary, expandable to the power series function f(x) = 
S^Lo a nX n i then the formal series 

oo 

f(x,t) = e» tD *f(x) = ]Ta„e^x" (24) 

OO 

= y j a n H N (x,vt;q), (25) 



represents a time dependent solution of the q-heat equation ()13[) The domain 
of convergency for this series is determined by asymptotic properties of our 
q-Kampe-de Feriet polynomials for n — > oo and requires additional study. 

According to this we have the evolution operator for the q-heat equation as 

U(t) = e vtD *. (26) 
It allows us to solve the initial value problem 

(^- V Dl)4>{x,t) = 0, (27) 
0(o;, 0) = f{x), (28) 

in the form 

ct>(x,t)=e» w *<t>(x,0)=e» tD lf(x), (29) 
where we imply the base q > 1 so that e q (x) is an entire function. 
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6 q-Burgers' Type Equation 

By using the q-Cole-Hopf transformation [9] 

u(x,t)=-2v , (30) 

(p(x,t) 

where <fi(x,t) is a solution of the q-heat equation (|13p. we find then that u{x,t) 
satisfies the q-Burgers' type evolution equation with cubic nonlinearity 

§- t u{x,t) - uD 2 x u{x,t) = i [(1 - M*)u(x,t)D x u{x,t)] - 

| [D x (u(qx, t)u(x, t))] + ± [u(q 2 x, t) - u(x, qt)] u(qx, t)u(x, t), 

where is the delation operator M*f(x) = f(qx). When q 1 this equation 
reduces to the standars Burgers' Equation 

u t + uu x = vu xx (31) 
6.1 I.V.P. for q-Burgers' Type Equation 

Substituting the operator solution (|29|) to (|3Q|l we find operator solution for the 
q-Burgers type equation in the form 



«(x,t) = -2 ^„^r;;_\ - ■ (32) 



^D x f{x ) 
e vtD lf{x) 

This solution corresponds to the initial function 



u(x,0) = -2vM^. (33) 
fix) 

Thus, for arbitrary initial value u(x,0) — F(x) for the q-Burgers equation we 
need to solve the initial value problem for the q-heat equation (|13[) with initial 
function f(x) satisfying the first order q-difference equation 

(D x + ^F(x))f(x) = 0. (34) 

7 q-Shock soliton solution 

As a particular solution of the q-heat equation we choose first 

4>(x,t) = e k2t e q (kx) , (35) 
then we find solution of the q-Burgers equation as a constant 

u(x,t) = -2vk. (36) 
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We notice that for this solution of the q-heat equation, we have an infinite set 
of zeros, and the space position of zeros is fixed during time evolution at points 
x n = -q n+1 /(q-l)k, n = 0,l,.... 

If we choose the linear superposition 

<f>(x, t) = e fc i*e g (kix) + e*»*e, (k 2 x) , (37) 
then we have the q-Shock soliton solution 

, . _ n he^eg {k lX ) + k 2 e k ^e q jk 2 x) 

u(x,t) — -Zv p- — — p- — — . (66) 

e K i l e q ikix) + e 2 e q ik 2 x) 

This expression is the q-analog of the Burgers shock soliton and for q — > 1 it 
reduces to the last one. However, in contrast to the standard Burgers case, 
due to zeroes of the q-exponential function this expression admits singularities 
coming from x for some values of parameters k\ and A^.To have regular solution 
we can follow similar approach as discussed in [S], then for k 2 = —k\ we have 
the stationary shock soliton 

, Ga(kix) — e a (—kix) „ , sinfuffciir) „ , , . 
e q \k\x) + e q (— kix) cosh g (fcix) 

(39) 

This function has no singularities on the real axis x and everywhere we have 
regular q-shock soliton as in [5] .However time evolution of shock solitons in [5] 
produce singularity at finite time. Here we like to find regular in x shock soliton 
which is regular at any time. 

We can choose solution of q-Heat equation (fT3")) as 

<p{x, t) = 10 + e fc i*e, (k x x) + e k ^e q (k 2 x) , 
then for k\ = 1 and k 2 — — 1 we get the q-Shock soliton 

e q (x) - e q (—x) 



u(x, t) = —1v 



lOe * + e q (x) + e q {— x) 



This solution describes evolution of shock soliton, so that at t —> — oo, u(x, t) — > 
0, and for t — > oo, u(x,t) — > — 2i/tanh 9 x In Figures 1,2,3 we plot the regular 
q-shock soliton for k\ = 1 and k 2 = —1 at different time t = —2, 0, 5 with base 
q = 10. 

If we plot the regular q-shock soliton evolution for k\ = 1 and k 2 = — 1 at 
different ranges of x and with q = 10, it is remarkable fact that the structure of 
our q-shock soliton shows self-similar property in the space coordinate x. Indeed 
at the range of parameter —50 < x < 50, and —5000 < x < 5000, structure of 
shock looks almost the same. 

For the set of arbitrary numbers ki, fc/v 

N 

(t){x,t) = y £ j e k ^ t e q {k n x), (40) 

71=1 
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Out[3]= 
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Figure 1: q-shock evolution for v = 1, k\ = 1, &2 = —1, t — —2 at range (-50, 
50) 



we have multi-shock solution in the form 



u(x, t) = —2v 



V w k 



€ n €q [k n X] 



iV 



e fe »*e 9 (k n x) 



(41) 



In general this solution admits several singularities. To have a regular multi- 
shock solution we can consider the even number of terms N = 2k with opposite 
wave numbers. When N = 4 and ki = 1, &2 = —1^3 = 2,^4 = —2 we have 
q-multi-shock soliton solution, 



u(x, t) = —2v 



sinh 9 (a;) + 2e A sinh 9 (2x) 
cosh g (x) + e 3 * cosh g (2x) 



(42) 



In Figures 7,8,9 we plot N — 4 case with values of the wave numbers k\ = 1, 
&2 = —1, &3 = 2, &4 = —2 at i = —10,0, 7 and with q = 10. This multi-shock 
soliton is regular everywhere in x for arbitrary time t. This result takes place 
due to absence of zeros for the standard exponential function e fc *. 

If we plot this regular multi q-shock soliton evolution at different ranges of 
x, it is remarkable that the structure of this regular multi q-shock soliton shows 
also the self-similar property in the space coordinate x. Indeed at the range of 
parameter —50 < x < 50, and —5000 < x < 5000, structure of multi q-shock 
looks almost the same. 
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Figure 2: q-shock evolution for v = 1, kx = 1, &2 = — 1, f = Q at range (-50, 50) 

8 Standard Time-Dependent q-Schrodinger Equa- 
tion 

The above consideration can be extended to the time dependent Schodingcr 
equation with q-dcformcd dispersion. We consider the standard time-dependent 
q-Schodinger equation 

!-»!*) (43) 

where ip(x,t) is a complex wave function. 
One can easily see that 



tp(x,t) = e %h2mt e q [-px 



i 

n 1 

is the plane wave solution of (I43|) . By expanding this in terms of momentum p 

'Tl' \N} : 

we get the set of complex q-Kampe-de Feriet polynomial solutions 



H { °\x,it;q) = Y, 
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Figure 3: q-shock evolution for v = 1, ki = 1, k% = —1, t = 5 at range (-50, 50) 
for O. 

Let us consider the complex version of the q-Cole-Hopf transformation 



u(x, t) = — 



ih D x ip(x, t) 
m if){x,t) 



then complex velocity function u(x, t) satisfies the complex q-Burgers-Madelung 
type equation 

ih§- t u(x,t) + ^D 2 x u{x,t) = fu(x,t)[l - M£]D x u(x,t)- 



ih 



[D x (u(qx, t)u(x, t))] + y [u(q 2 x, t) — u(x, t)~\ u(qx, t)u(x, t). 



If we write u — u± + iu 2 and separate it into real and imaginary parts , then 
we get two fluid model representation where u\ is a Madelung-London-Landau 
classical velocity, and u 2 is the quantum velocity. 
For the real part we have 

d ft 2 if i 

~h-7^u 2 {x,t) + —D 2 x ui(x,t) = — [(u^qx^) - ui(x, (x,t)ui (qx,t) - 



u 2 (x,t)u 2 (qx,t)) - (u 2 (q 2 x, t) - u 2 (x, t))(m (x, t)u 2 (qx, t) + u 2 (x, t)m(qx, t))) 
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Figure 4: q-multi shock evolution for ki = 1, k 2 
and at range (-50, 50) 



■l,k 3 = 2,k 4 = -2, t = -10 



-~[ui(s,t)[l - Af?]D s « 2 (»,t) +ua(x,t)[l - M?\D x ui(x,t)\ 



+ -D x [u 2 (qx,t)u 1 (x,t) + ux(qx,t)u 2 {x,t)], 



(44) 



and for the imaginary part 

d fv TTl 

h—uxix.t) + —D 2 x u 2 {x 1 t) = —[(u 1 (q 2 x,t) - ui{x,t)){ui{x,t)u 2 {qx,t) 
at 2m 2 

+u 2 (x,t)ui(qx,t)) + (u 2 (q 2 x 1 t) - u 2 (x,t))(u\(x,t)u\(qx,t) - u 2 (x, t)u 2 (qx, t))] + 



-[ Ul {x,t)[l - M*]D xUl (x,t) - u 2 (x,t)[l - M x q ]D x u 2 {x,t)\ 



-D x [ui(qx,t)u 1 (x,t) - u 2 (qx,t)u 2 {x,t)]. 
When q — > 1, the real part reduces to the continuity equation 



(45) 



-{u 2 ) t + —{ui) xx = (uiW 2 )x, 
2m 
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Figure 5: q-multi shock evolution for ki — 1, &2 — —1,^3 = 2, k^ = —2, t = 
and at range (-50, 50) 

and the imaginary part reduces to the Quantum Hamilton- Jacobi equation 



Thus the two fluid system (|4"4")l . (|45p is the q-analogue of the coupled q- 
quantum Hamilton- Jacobi equation and the q-continuity equation. 

Following similar procedure as in first part of this paper, we can construct 
particular solutions of our q-Schrodinger equation in the form of complex shock 
solitons. This question is under investigation now. 
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